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Abstract. Single spin asymmetries in semi-inclusive deep-inelastic scattering off transversely po- 
larized targets give information on, among other fascinating effects, a pseudo time-reversal odd 
parton distribution function, the 'Sivers function'. In this proceeding 1 we review the extractions of 
this function from HERMES and COMPASS data. In particular, the HERMES pion and kaon data 
suggest significant sea-quarks contributions at x ~ 0. 15 to the Sivers effect. We present a new fit that 
includes all relevant sea quark distributions and gives a statistically satisfactory overall description 
of the data, but does not describe ideally the K + data from HERMES. We argue that measurements 
of the pion- and kaon Sivers effect at CLAS12, and COMPASS, will clarify the situation. 
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1. INTRODUCTION 

Ever since the first large single spin asymmetries (SSA) were observed in hadron-hadron 
collisions [1, 2, 3, 4] spin phenomena in QCD became more and more popular. In par- 
ticular the SSA in semi-inclusive deeply inelastic scattering (SIDIS) from transversely 
polarized targets have recently been measured at HERMES and COMPASS [5, 6, 7, 8, 
9, 10, 11, 12, 13, 14]. SSA in SIDIS with longitudinally polarized leptons or nucleons 
have been reported at HERMES and the Jefferson Lab [15, 16, 17, 18, 19, 20, 21]. In a 
partonic picture structure functions in SIDIS and the Drell-Yan (DY) process with small 
transverse momenta P± <C Q of the final state hadron or the lepton pair, respectively, 
are sensitive to the intrinsic transverse motion of partons [22, 23, 24], see also the re- 
view [25]. This means in the case of SIDIS that structure functions can be expressed in 
terms of convolutions of transverse momentum dependent (TMD) parton distributions 
and fragmentation functions. Strict factorization formulae for these processes involving 
an additional soft factor were discussed and established in [26, 27, 28]. 

In the production of unpolarized hadrons in DIS of an unpolarized lepton beam off 
a transversely polarized target there are three leading twist spin structure functions, see 
e.g. [25], which can be distinguished by their azimuthal distributions proportional to 
sin(0 + S ), sin(0 — S ), and sin(30 — S ). Here denotes the angle between lepton 
plane and hadron plane in the lab frame while S represents the angle of the transverse 
spin vector with respect to the lepton plane. In the parton model the structure function 
proportional to sin(0 + $$) ("Collins effect") is expressed in terms of two chirally-odd 
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FIGURE 1. The SSA due to Sivers effect arises in SIDIS from final state interactions [33] (upper 
part), and in DY from initial state interactions [34] (lower part). Both types of interactions are encoded 
appropriately defined Wilson lines that are connected to each other by time reversal [35] . In the Figure the 
respective interactions are sketched in the one-gluon-exchange approximation, see text. 

correlation functions, namely the transversity parton distribution h\ which is not acces- 
sible in inclusive DIS, and the Collins fragmentation function [29]. The structure 
function proportional to sin(0 — ty s ) ("Sivers effect") is described by the T-odd parton 
distribution /[^, the so-called Sivers function [30, 31], in conjunction with the usual un- 
polarized fragmentation function D\ . The structure function proportional to sin(30 — (j) s ) 
provides information on the so-called "pretzelosity" distribution h\ T [23] whose physical 
interpretation was recently discussed in [32]. 

T-odd parton distributions such as the Sivers and also Boer-Mulders function h\ [24] 
were considered to vanish due to time-reversal symmetry for some time. Only when 
it became clear that initial/final state interactions between the struck quark and the 
target remnants in the parton model can cause SSAs [33, 34], see Fig. 1, the existing 
definitions of transverse momentum dependent parton distributions were revisited. It was 
then shown that initial/final state interactions can be implemented into the definitions 
of parton distributions by means of gauge-link operators with appropriate Wilson lines 
[35, 36, 37, 38]. In particular these Wilson lines ensure the color gauge invariance of 
the definition of TMD parton distributions and fragmentation functions. The application 
of time-reversal switches the direction of the Wilson lines from future-pointing (final 
state interactions) to past-pointing (initial state interactions) lines in conjunction with an 
overall sign change. Therefore T-odd parton distributions do not vanish. Instead, time- 
reversal establishes a connection between T-odd parton distributions in processes with 
final state interaction (e.g. SIDIS) and in processes with initial state interactions (e.g. 
Drell-Yan). For the Sivers function the relation reads 



This important QCD-prediction can and still needs to be checked by experiments. 




SIDIS 




DY 



(1) 



In these proceedings we discuss the present status of the understanding of the Sivers 
function from HERMES and COMPASS [39, 40, 41, 42, 43, 44, 45, 46, 47]. In particular 
we review the works [39, 42, 43, 44, 46, 47] and discuss the developments due to the 
most recent data from SIDIS. We also present predictions for pion and kaon Sivers 
asymmetries for the 12 GeV upgrade at Jefferson Lab. 



2. SEMI-INCLUSIVE PROCESSES 

In the following we briefly discuss SIDIS with transversely polarized targets, with 
particular emphasis on the Sivers effect. In SIDIS of an unpolarized lepton beam off 
a transversely polarized nucleon target one can measure the SSAs [48, 49, 50, 25] 
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The SSAs A^p^ = F^f^/Fuu are ratios of the respective structure functions to the 

unpolarized one, which depend on x = Q 2 / (IP ■ q) , z = (P ■ Ph) / (-P • <?) > P%± an ^ Q 2 = ~Q 2 
where P is the target momentum. The angles and other momenta are defined in the left 
panel of Fig. 2, and £ ~ (1 — y)/(l —y + y 2 /2) denotes the polarization of the virtual 
photon with y = (P ■ q) /(P ■ I), and St is the transverse spin vector. The dots in Eq. (2) 
denote terms with other angular distributions due to pretzelosity or subleading twist. 

The parton model expressions for the structure functions in terms of convolutions of 
TMD parton distributions and fragmentation functions can be easily understood when 
using the "tree-level" formalism of Ref. [23] where only the leading order (in oc s ) tree 
diagrams in the hard part are considered. By splitting the SIDIS cross section into 
a leptonic tensor L^ v (assuming one-photon exchange) and a hadronic tensor W^ v , 
<^Osidis 06 L^ V W^ V , we can express the leading part (in 1/0 of W^ v by the tree- 
diagram in the right panel Fig. 2. The lower soft blob in that diagram describes the 
momentum distribution of quarks inside the target whereas the upper part represents the 
fragmentation of a parton into hadrons. Field theoretically the soft blobs correspond to 




FIGURE 2. Left panel: The kinematics of the SIDIS process IN — > I'hX. The nucleon is polarized 
transversely with respect to the beam. However, up to power corrections the polarization is transverse also 
with respect to the momentum of the virtual photon: sin ©5 ~ M;v /Q <C 1 . Right panel: Tree-level diagram 
for the hadronic tensor in leading order in 1 /Q in the parton model. 



the following matrix elements of bilocal quark-quark operators, 
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The gauge link operators and m Eqs. (3) and (4) ensure the color gauge invari- 
ance of the matrix elements. They are defined as 

Wct>[0, (z~ ,0,zt)\v] = [0|a*] x [a$,\b<t>] x | c<*>] x [c<t| (z~,0,zr)], (5) 
^ A [oo,(0,z+,z r )|77] = [& A |c A ]x[c A |(0,z+ z r )], (6) 

#A[0,°o|Tf] = [0|flA]x[fl A |fcA], (V) 

where [alb] denotes a gauge link operator, i.e. a path-ordered exponential of gluon 
field operators, with a straight Wilson line between the space-time coordinates a and 
b. The three "milestones" of the Wilson lines in Eqs. (5)-(7) are a<j> = (r]oo,0,07-), b$> — 
(r]oo 5 0,°<>r) and c<p = (T7°o,0,?r) for the correlator 4>, while we have a A = (0, r/°o,0r), 

= (0, T7°°,°°r) and c A = (0, r/°o 7 J r ) for the fragmentation correlator A. Therein we 
use the usual light cone coordinates for a four vector a} 1 = (a^,a + ,ar) with = 

(a° ± a 3 ) . The parameter r\ e { — 1,1} determines the direction of the Wilson line. As 
was pointed out in [35] time-reversal transforms a Wilson line with r\ = +1 describing 
final state interactions, a situation one encounters in SIDIS, into a Wilson line with 
r\ = — 1 describing initial state interaction, e.g. in Drell-Yan. This results in the relation 
(1) for the Sivers function. 

It is convenient to express the diagram in the right panel of Fig. 2 in terms of the 
correlators (3) and (4) in a frame in which the momentum of the nucleon P and the 
momentum of the produced hadron P^ are collinear. We obtain for the leading part in 
1 /Q of the hadronic tensor 
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The quark-quark correlators 4> and A contain all possible spin information of either the 
probed or fragmenting quark as well as the spin information of the target. It is possible 
to project out various polarizations of the quarks by tracing the correlators with appro- 
priate Dirac matrices. These traces can be then parameterized in terms of TMD parton 
distributions and fragmentation functions [23, 24] (for a complete parameterization see 
[51]). Assuming that the nucleon moves fast in one light cone direction, i.e. P + is the 
large component of the nucleon momentum, we obtain two TMD distributions for un- 
polarized quarks by tracing 4> with y + . One is the well-known ordinary PDF for an 



unpolarized target f\ while the other is the Sivers function, describing the distribution of 
unpolarized quarks in a transversely polarized nucleon (e^ = £~ + ' J ' with £ 0123 = +1) 

^Tr[<J>(*,p r )y+] = fl (x,p 2 T ) - ^^ftAxS)- (9) 

By tracing 4> with other Dirac structures we obtain further PDFs. For example, by tracing 
it with i(7 l+ f5 we obtain, among others, the transversity distribution hi and the so-called 
'pretzelosity' distribution hf T , etc. 

Since here we consider unpolarized hadrons in the final state only, there are only two 
relevant fragmentation functions, the ordinary fragmentation of unpolarized quarks and 
the chirally-odd Collins function of transversely polarized quarks. Assuming that the 
produced hadron moves fast in the minus light cone direction, i.e. the large component 
of Ph is P^, we project them out by tracing with y~ and io l ^ Ys , 

1 1 E^k? 

-Tr[A( z ,* r )y-] = Dfaft) ; -Tr[A(z,^)io , - 1S ] = — ^^(z,^). (10) 

At this point one obtains the expression for the structure functions Sivers structure 
function by inserting the traces of the type (9), (10) into the hadronic tensor (8), boosting 
into a frame in which the momentum of the nucleon and of the virtual photon are 
collinear, and contracting with the leptonic tensor L^ v which yields: 
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where h = Ph±/\Ph±\ I s the normalized transverse momentum of the produced hadron. 
The unpolarized structure function Fyy is given by (1 1) with (—h-p T ) f^- q /M h — > ff. 

We emphasize again that the result (11) is valid for the tree-level hard process. A 
factorization valid to all orders requires the inclusion of a soft factor [26, 27, 28] in the 
formula (11) in order to handle soft gluon radiation. Up to now no phenomenological 
treatment takes this factor into account. We will neglect this factor in the following. 



3. SSA IN p T p -> ttX 



The first information on the Sivers function was obtained from studies [52, 53] of SSAs 
in pip — > nX or jfip — > nX [3, 4]. Although they originally motivated the introduction 
of the Sivers effect [30, 31] the theoretical understanding of these processes is more 
involved compared to SIDIS or DY. Here SSAs can also be generated by twist- 3 effects 
[54, 55, 56, 57, 58], though it was suggested that these could be manifestations of 
the same effect in different kj regions [59, 60, 61]. Moreover, for this processes no 
factorization proof is formulated in terms of the Sivers effect. Studies of other processes 
with pp in the initial state and hadronic final states indicate that it is in general hard to 
prove factorization in hadronic reactions [59, 62, 63]. 



4. SIVERS EFFECT IN SIDIS: FIRST INSIGHTS 



The first information on the Sivers function from SIDIS was obtained in [39] from a 
study of preliminary HERMES data [64] on the 'weighted' SSA defined as 
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where N- are sums over event counts for the respective transverse target polarization, 
and (...) denotes averaging — here over z and The advantage of 'weighted SSAs' 
is that the integrals in the structure function (11) can be solved exactly [24] yielding 

A P h±/ M NM ^s) ( , _ 2fdP^F^\ X ^) _ (-2) Laelxf^WDUz) 
l ' J ' fdP^Fuvix^) ' L a elxf?(x)D«(z) 

(13) 

where f^ a (x) = Jd 2 p T jjfrfjf(x,ft). 

While the weighting is preferable from a theory point of view, it makes data analysis 
harder. It is difficult to control acceptance effects, and the HERMES Collaboration does 
not recommend the use of the preliminary data [64]. In 'unweighted SSAs' defined as 
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acceptance effects largely cancel. Therefore such data have been finalized first, and one 
even is not discouraged to use preliminary data of this type [7, 9, 10]. However, the prize 
to pay is that now the convolution integrals in (11) can be solved only by resorting to 
models for the transverse momentum dependence. Here we assume the distributions of 
transverse parton and hadron momenta in distribution and fragmentation functions to be 
Gaussian with the corresponding Gaussian widths, p| iv and K% , taken to be x- or z- and 
flavor-independent. The Sivers SSA (14) as measured in [5, 6] is then given by [42] 
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In view of the sizeable error bars of the first data it was necessary to minimize the number 
of fit parameters. For that in [42] effects of sea quarks were neglected. In addition, the 
prediction from the limit of a large number of colors N c in QCD [65], namely 

frr(x,pj-) = -fij{x : p T ) modulo \/N c corrections, (16) 

was imposed as an exact constraint. Analog relations holds also for antiquarks, and all 
are valid for x of the order xN c = 0{N^) [65]. The following Ansatz was made and 
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FIGURE 3. (a) Sivers function of M-quarks vs. x at a scale of 2.5 GeV, as obtained from HERMES 
data [5]. Shown are the best fit and its 1-c uncertainty, (b + c) The Sivers SSA, Eqs. (14, 15), for n + from 
proton as function of x and z as obtained from the fit in Figure 3a in comparison to the data [5]. The n~ 
data from [5] are compatible with zero, and are equally well described (not shown here). 



best fit obtained: xf^ X)u {x) = -xf^ X) \x) A = tz Ax b (1 -x) 5 = -0.17x a66 (l -x) 5 
[42, 43]. Fig. 3a shows the fit and its 1-a uncertainty due to the statistical error of the 
data [5]. Fig. 3b shows that this fit well describes the x-dependence of the HERMES data 
[5]. Fig. 3c finally shows the equally good description of the z-dependence of the data 
[5], that were not included in the fit. This serves as a cross check for the Gauss Ansatz, 
which apparently works well here. It is found in general that this model is useful as long 
as one deals with limited precision and small transverse momenta (Phi) <C Q [66, 53]. 

In [42] it was furthermore found that effects due to sea quarks could not be resolved 
within the error bars of the data [5]. It was also checked that 1 /A^ c -corrections expected 
in (16) are within the error bars of HERMES [5] and especially COMPASS [8] data. 

To draw an intermediate conclusion, the HERMES and COMPASS data [5, 6] are 
compatible with large-/V c predictions [65] . More precisely, the large-/V c approach worked 
at that stage, because the precision of the first data [5, 6] was comparable to the theoret- 
ical accuracy of the large-A^- relation (16). Remarkably, the signs of the extracted Sivers 
functions, fyp < and fyf > 0, agree with the physical picture discussed in [67]. The 
findings of [42] were in agreement with other studies [40, 41], see also the review [45]. 



5. SIVERS EFFECT IN DY: FIRST PREDICTIONS 

As the experimental test of the particular 'universality relation' for the Sivers function in 
Eq. (1) is of fundamental importance, the first insights on f^. [39, 42] were immediately 
used to estimate the feasibility of DY-experiments to measure the Sivers effect [39, 42]. 
It was shown that the Sivers effect leads sizeable SSAs in p^n~ — > l + l~X, which can 
be measured at COMPASS, and in p or pp^ — > l + l~X, which could be studied in the 
proposed PAX experiment at FAIR, GSI [68]. 

On a shorter term the Sivers effect in DY can be studied in p^ p — > l + l X at RHIC. 
In pp-collisions sea quarks are involved, and counting rates are smaller. It was shown, 
however, that the Sivers SSA in DY can nevertheless be measured at RHIC with an 
accuracy sufficient to unambiguously test Eq. (1) for Sivers quarks [42, 44]. Moreover, 
RHIC can provide valuable information on Sivers antiquarks [42, 44]. 



6. SIVERS EFFECT IN SIDIS: FURTHER DEVELOPMENTS 



The first data [5, 6] gave rise to a certain optimism. The Sivers effect on a proton target 
was clearly seen [5], its smallness on a deuteron target [6] understood in theory [65]. 
Independent studies [39, 40, 41, 42] agreed on the interpretation of the SIDIS data [45]. 
Prospects of accessing the effect in DY (and to test the universality property of f^f) 
were found promising, see [39, 40, 41, 43] and [69, 70, 71]. The optimism persisted the 
following releases of higher statistics data [7, 8] from HERMES and COMPASS, but it 
was somehow damped with the advent of HERMES data on kaon Sivers SSAs [9, 10], 
while COMPASS reconfirmed that on deuteron the effect is really small [11, 12, 13]. 

Let us first discuss what one would naively expect for the K + Sivers effect. Comparing 
the 'valence quark structure' (K + ~ us vs. 7t + ~ ud), we see that the DIS -production of 
K + and n + differs by the 'exchange' of the sea quarks s <-> d. If it were legitimate to 
neglect sea quarks in nucleon, then the K + and n + Sivers SSAs would be comparable. 

Now let us have a look on the preliminary HERMES data [9]. At larger x > 0.15 
the Sivers SSAs for K + and n + are comparable, as expected. The fit to pion data 
discussed in Sec. 4 well describes the K + data in this region [46, 47], see the solid line 
in Fig. 4a (where, of course, all 'differences' between kaons and pions due to different 
fragmentation functions [72] are considered). However, at smaller x ~ 0.1 we observe a 
(2-3) times larger SSA for K + compared to n + . For K ~ su (i.e. pure 'sea quark effect') 
the SSA is compatible with zero and bears no surprises. 

The K + Sivers effect at HERMES [9, 10] hints at an importance of Sivers sea quarks. 
Interestingly, no kaon over pion enhancement is observed at COMPASS [11, 12, 13]. 
This is remarkable, because the COMPASS kinematics covers much smaller x down to 
x m i n = 0.003 (vs. HERMES x m [ n = 0.023). Sea quark effects could therefore show up at 
COMPASS even more clearly — however, one has to keep in mind the different targets: 
proton at HERMES vs. deuteron at COMPASS. 

In order to get a feeling about the impact of sea quark effects, let us do the following 
exercise. We use for f^f and fyf the previous best fit results, see Sec. 4, and add on top 
of that Sivers u, d, s and ^-distributions which saturate the positivity bounds [73] 

l/f°WI<^W, (iv) 

where (pf) U np is the mean transverse momentum of unpolarized quarks in the nucleon 
[42]. (It could depend on flavor. Neglecting this possibility one obtains a good des- 
cription of transverse hadron momenta at HERMES [17] for (pr) un p = 0.5 GeV [42].) 
The effects of sea quarks allowed to saturate ± these bounds are shown in Fig. 4a as 
dashed lines. Throughout the parameterization [74] is used for f". 

Fig. 4a demonstrates that sea quarks may have a strong impact, and could be able 
to explain the Sivers K + SSA. However, introducing Sivers sea quarks as large as to 
explain the K + data, Fig. 4a, tends to overshoot the n + data, see Fig. 4b. Also there is 
no reason to expect Sivers sea quarks to saturate positivity bounds. 

If the Sivers K + effect is not a statistical fluctuation (which is unlikely, see Fig. 4c), 
then it should be possible to obtain a satisfactory fit to the pion and kaon data from 
HERMES and COMPASS. The next Section is devoted to this task. 
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FIGURE 4. (a) The Sivers SSA for K + as function of x. The preliminary HERMES data is from [9] . The 
solid line is the K + SSA obtained from the best fit to pion data [5] (see Sec. 4). The dashed lines display 
the effect of adding on top of that Sivers sea quarks saturating ± the positivity bounds (see Sec. 6). It seems 
that sea quarks could explain the effect, but at the same time one overshoots n + data (see next figure), 
(b) Sivers SSA for n + as function of x. The published HERMES data are from [5], and the theoretical 
curves as in Fig. 4b (but for n + ). (c) Comparison of the first (lower statistics, boxes) [9] and the most 
recent (higher statistics, circles) [10] data from HERMES on the K + Sivers SSA as function of x. It seems 
unlikely that the effect could be due to statistical fluctuation, especially in the region of x ~ 0.1. 



7. UNDERSTANDING PION AND KAON SIVERS EFFECT 

In order see whether it is possible to understand the data on the Sivers SSAs for pions 
and kaons from different targets, it is necessary to attempt a simultaneous fit. We use 
the HERMES proton target data on 7t°, 7t ± , if ± from [10] and the COMPASS deuteron 
data on n ± , from [11]. Since we do not know error correlation matrices we cannot 
perform a simultaneous fit to data on x, z and P/ 7 ± -dependences with correct estimate of 
its statistical significance. In the present study data on jc-dependence will serve as input 
for the fit, and data on z-dependence will be used only for a cross check of the results. 
When using data from different experiments it is necessary to consider systematic 
uncertainties. Those are dominated by the uncertainty of the target polarization in both 
experiments, and we combine them with statistical errors in quadrature. 

The data on x-dependences of the Sivers SSAs have little sensitivity to the parameters 
entering the factor ao in (15). We fix the Gaussian width of the fragmentation function 
and Sivers function as (Kjj ) = 0.16GeV 2 [42], and (p| iv ) = 0.2 GeV 2 which coincides 
with the central value obtained in [42]. We stress that the final results presented here 
anyway depend only weakly on these parameters, which could be inferred from studies 
of data on /^-dependence. Such studies will be reported elsewhere. Finally, for f"(x) 
and D"(z) we use the leading order parameterizations [72, 74] at a scale of 2.5 GeV 2 
which is close to the (Q 2 ) of both experiments. 

Now we need an Ansatz for the Sivers functions. Here we shall content ourselves to 
the following simple Ansatz (the constraints on the A a arise from positivity, Eq. (17)) 

fir 1)a =^ { -^ff(x), \A a \<l (18) 

with (pr)unp = 0.5 GeV from [42] . Since sea quark effects are of importance (see Sec. 6) 
we introduce the flavours: a = u,d, s, w, J, s. So the initial task is to fix the A^ pa ra = 6 
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FIGURE 5. The Sivers SSA for various hadrons from different targets vs. x. Left panel: HERMES data 
(proton target) [10]. Right panel: COMPASS data (deuteron target) [11]. The theoretical curves are the best 
fit (solid line) and its 1-c-region (shaded area) as obtained from the Ansatz and best fit in Eqs. (18, 19). 
These data served as INPUT for the fit (18, 19). 



parameters A a from A^ata = Mermes + A/compass = 25 + 36 data points. The Ansatz 
(18) makes the numerical handling of the problem particularly simple. The x 2 is a 'six- 
dimensional parabola' in the space of the A a , and the only extremum (global minimum) 
is easily found by means of simple self-made codes, or minuit [75]. 

The best fit has a x 2 = 73. 4. This means a satisfactory x 2 P er degree of freedom of 
Z 2 /C^data — A/para) = X 2 f = 1-33. The results for the best fit parameters read: A u = 
-0.21, A d = 0.41, A a = 0.24, A d = -0.27, A s = 0.95, A s - = -1.93. The (correlated!) 
1-a uncertainties of these fit results are of ^(10%). But A $ exceeds the positivity bound 
(18) and A s comes suspiciously close to it, which is driven by the K + HERMES data. 

We therefore repeat the fit and use the above results to inspire the following Ansatz. 
We fix A s = +1 and A.y = —1 from the very beginning. This 4-parameter-fit has, of 
course, a slightly higher x 2 = 76.5 but nearly the same Xdof = 1-34 which means that 
it is equally good. The best fit parameters read 

A u = -0.21 ±0.01 A a = 0.23 ±0.02 A, fi = d ±l 

A d = 0.38±0.03 A d = -0.28 ±0.04 A. f fi = d -1. (19) 

The 1-c errors are those solutions for the A a (a = u, d, u, d) which increase the total 
X 2 = 76.5 by one unit. For our Ansatz this is again particularly simple. For example, the 
uncertainty of A u is found by fixing the other free parameters to their best fit values, and 
solving a quadratic equation. The 1-a uncertainties in (19) are correlated, of course. 

Fig. 5 shows how the fit (18, 19) describes the data on the ^-dependence of the Sivers 
SSAs (i.e. the input for the fit). Let us comment on the fit: 

• The fit quality is satisfactory: X 2 f = 1-3 = ^(1) as it should be. (There are 
stronger criteria for goodness of a fit [76], but we have too few data to apply them.) 

• Looking at Fig. 5 we notice: out of 61 data points, there are only two(!) data points 
that are off the best fit curve in a worthwhile mentioning way. 

• From the point of view of statistics, one can comfortably live with such a situation, 
and wait for new data that will allow to improve the fit. 

• But is it not suspicious that those two points, that are off the best fit in a worthwhile 
mentioning way, are precisely the K + HERMES data around x ~ 0.1? 

The last observation raises the question, whether the K + HERMES data around x ~ 0. 1 
could be a statistical fluctuation. Of course, this possibility cannot be excluded — though 
one is not inspired to consider such an explanation as convincing, looking back at Fig. 4c. 

Let us postpone the discussion of this point for later, and continue with z-dependence. 
We did not use these data in the fit. Therefore the results, see Fig. 6, are a prediction of 
the fit (18, 19), and the assumed Gaussian model for transverse parton momenta. The 
overall performance of the description is again satisfactory. It is not surprising that the 
only worthwhile mentioning mismatch is for the K + data from HERMES. 

At this point, one could try to improve the description of z-dependence by re-adjusting 
the parameters and in Eq. (15). In particular, could even be allowed to be z- 

dependent. With different and Kp one would, of course, obtain somehow different 
A a in the best fit (18, 19). This could be continued to an iteration procedure, which (if 
convergent) would result in a optimized description of data on x- and z-dependence. 



. sin(4"-4> s > 
UT 



HERMES preliminary 



COMPASS preliminary 



0.15 - 

0.1 - 

0.05 r 

r 

-0.05 - 

-0.1 - 



.sln(i|H|i s ) 
UT 



prolon 7t 



vs. HERMES preliminary 



0.15 - 

0.1 - 

0.05 } 

r 

-0.05 - 

-0.1 — 



A UT 



prolon Tl 



~1 



HERMES preliminary 



0.15 - 
0.1 

0.05 

- 

-0.05 ~ 7 

-0.1 - 



prolon K + 





A UT 



0.15 
0.1 
0.05 



HERMES preliminary 



A UT 



COMPASS preliminary 



-0.05 
-0.1 



UT 



0.15 
0.1 
0.05 



i .s. HERMES preliminary 



: proton w 






L-i- \ 


- , , , , i , , , , i , , , , i , , , , - 



_ i i i i | i i i i | i i i i | i i i i 
: proton K~ 


0.1 
0.05 


_ i i i i | i i i i | i i i i | i i i i | i 

: deuleron K~ 















: t-riirL 






-0.05 








-0.1 






- 


-0.15 


- , , . , i , . . , i , . , , i . , 





0.2 0.4 0.6 0.8 



FIGURE 6. The Sivers SSA for various hadrons from different targets vs. z. Left panel: HERMES data 
(proton target) [10]. Right panel: COMPASS data (deuteron target) [11]. The theoretical curves are the best 
fit (solid line) and its 1-c-region (shaded area) as obtained from the Ansatz and best fit in Eqs. (18, 19). 
These results are a PREDICTION of the fit (18, 19). 



Phi. -dependence could also be included. This procedure would help to better constrain 
the parameters p^ iy and K^^ . But it presumably would have little impact on the fit results 
(19). Keeping this in mind, we shall — at the present stage of art — content ourselves 
with the descriptions in Figs. 5 and 6. The Ph± -dependence will be discussed elsewhere. 

8. HOW DOES THE SIVERS FUNCTION LOOK LIKE? 

In Sec. 7 we have seen that the probably simplest Ansatz one can imagine for the Sivers 
function, see Eq. (18), works. 'It works' means in this context that it yields an acceptable 
X 2 per degree of freedom of Xd.o.f. = ^(1)- Other, more flexible Ansatze could yield 
better descriptions. In particular, it is an interesting question, whether one could better 
describe the proton target K + Sivers effect. Different Ansatze with more free parameters 
(and the impact of different fragmentation functions) were explored in [77], but a 'better' 
description of the kaon HERMES data could not be reported there, neither. 

However, having achieved a Xd.o.f. = 1-3 one must wonder, whether there really is a 
necessity to improve that fit. Let us adopt here the point of view that it is not, and draw 
conclusions from our results. 

• The Ansatz ( 1 8) is rather rigid. It denies the Sivers function an independent x-shape, 
and forces it to be proportional to /f (x). The only "freedom" it gives to f^f is that 
the proportionality factors A a can be flavor dependent. (To recall: we included the 
factor (pr)unp/ (2M#) in (18) such that \A a \ < 1 guarantees positivity (17).) 

• The initial, unconstrained six-parameter-fit forced the Sivers s (s) function to come 
close to (+1) x (to exceed ( — 1) x) the positivity bound (17). 

• Because of that, in the final four-parameter-fit, we fixed the Sivers strangeness 
functions such that they saturate ± the positivity bounds: A s = +1 and Ay = — 1. 

• One may worry whether such a large 'Sivers strangeness' in the nucleon is natural. 
However: (i) the 'net Sivers strangeness content' is zero (s, s have opposite signs), 
(ii) Fixing, for example, A s = j and A^ = — \ (explores positivity only within 50%), 
would increase the Xdof by only 0.07 units — i.e. an equally acceptable fit. 
Thus, presently there is no reason to worry about the Sivers strangeness functions. 
(Moreover (pr)unp in (17), defined within the Gauss model approximation, could 
be flavor-dependent and different for s, s, possibly relaxing numerically the bound 
imposed here.) 

• It is a reasonable guess that the Sivers function is suppressed at small x compared 
to fi(x) [41]. But it is not necessary to build in such a suppression in the Ansatz 
in order to describe the COMPASS data rather precise at small x. In our fit this is 
achieved by the different relative signs of the Sivers functions. 

• We do not mean that the Sivers function should raise like f^(x) at small x. Rather 
we would like to stress that a suppression of the Sivers function compared to f"(x), 
if existent, is not yet constrained by the present data. 

• We can draw from our study even the following stronger conclusion. The present 
data do not yet tell us much about the shape of the Sivers function. But they tell us 
something about its magnitude and relative signs of the different flavours. 



? l(l)a 




? l(l)a 



-0.02 



0.04 



0.02 




.l(l)a 



-0.02 



-0.04 
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After these cautious remarks concerning the meaning and interpretation of our results, 
let us have a look how the obtained Sivers functions look like, see Fig. 7. We make the 
following observations. 

• The Sivers u and d distributions are of comparable magnitudes but opposite signs, 
as predicted in the large-/V c limit [65], see Eq. (16). 

• The Sivers u and d distributions are of comparable magnitudes but opposite signs, 
as predicted in the large-/V c limit [65], see text in the sequence of Eq. (16). 

• It is 2A U ps — 2A a ~ A d ps — Aj and A s = —Ag, Eq. (19). Were the former relations 
exact, then the contribution of q, q (q = u, d,s) to the Burkardt sum rule [78, 79] 

£ fdxf$ l)a {x) = (20) 

would vanish, i.e. also Jq dx f^ 1 ^ 8 (x) = 0. That would imply a small Sivers gluon 
distribution, as concluded independently in [39, 80, 81]. 2 

• The Sivers u and d distributions obtained here agree qualitatively with earlier works 
[39, 40, 41, 42, 45] in which sea quarks effects were not considered. 

• The signs of the Sivers u and d distributions support the picture of Ref. [67]. 



2 Notice that strictly speaking the integrals in (20) exist only if A u = —A,-, and Aj = — Aj holds exactly 
in the Ansatz (18) (cf. in this context also the model calculation [82] where quark- and gluon Sivers 
functions cancel under the integral). For A u and Aq this is the case within 1 -o\ but for A c [ and A j there is 
a clearer mismatch, see (19), such that (20) would diverge. However, our parameterization is not meant to 
be valid for x — » and is anyway only constraint for x > 0.003 by COMPASS data. We remark that our 
parameterization also is not meant to be valid for jc — > 1 where the Sivers function is expected to vanish 
faster than f°(x) [83]. 



9. PREDICTIONS FOR CLAS AND COMPASS 



Fig. 8 shows the predictions, made on the basis of the fit (18, 19), for the x-dependence 
of the Sivers SSA of charged pions and kaons from proton targets for two experiments: 
CLAS with 12GeV beam upgrade, and COMPASS. The CLAS error projections for 
2000 hours run time are from [84, 85]. The COMPASS error projections are estimated 
assuming a statistics comparable to that of the deuteron target experiment [11]. 



shi((M>s) CLAS n GeV prediclion 



A 

0.15 - 



proton 71 



-0.05 - 

' ' 

0.1 0.2 0.3 0.4 0.5 0.6 * 



A UT Y Ts CLAS. 12 GeV, prediclion 



proton Tl 



0.1 0.2 0.3 0.4 0.5 0.6 x 



,sin(4>-(|> s ) 
UT 



COMPASS, prediction 



0.15 
0.1 



proton 1L 



T M i M } } 



COMPASS, prediction 




A UT Y Ts CLAS. 12 GeV, prediction 



0.1 0.2 0.3 0.4 0.5 0.6 x 



A UT 



COMPASS, prediclion 



0.15 
0.1 


T proton K + 


0.15 
0.1 


1 1 

T proton K + 




0.05 








0.05 





















iTTTl 







: i i [ [ 












0.05 






r , , , , 


-0.05 




■ . . i i 





UT 



0.15 

0.1 - 

0.05 - 

— 

-0.05 1 



CLAS, 12 GeV, prediction 



COMPASS, prediction 



proton 71 



0.1 0.2 0.3 0.4 0.5 0.6 x 




FIGURE 8. Predictions for the Sivers SSA as function of x from a proton target for CLAS at 12 GeV 
(left panel) and for COMPASS (right panel). The predictions are made on the basis of the fit (18, 19) 
presented in this work. The CLAS error projections are from [84, 85]. For the COMPASS error estimates 
we assumed the statistics of [1 1]. 



We see that both experiments will be able to confirm (or reject) these predictions 
within their expected statistical accuracies. 

The interesting question is, of course, whether CLAS and COMPASS will confirm 
the HERMES results for K + at x ~ 0. 15, and overshoot our predictions in that region of 
x. In this context, however, it is also worth to look more carefully into the kinematics. 
The next Section is devoted to this task. 



10. FACTORIZATION AND POWER CORRECTIONS 

When averaged over the respectively covered kinematical regions, CLAS, COMPASS 
and HERMES have a comparable (Q 2 ) ~ (2-3) GeV 2 (this is what we used in our study). 
At fixed x, however, Q 2 can vary significantly in these experiments. For example, 

HERMES: (x) =0.115, (Q 2 ) = 2.62 GeV 2 , 
COMPASS: (x) =0.1205, (Q 2 ) = 12.9 GeV 2 , (21) 

which is in the ^-region where the 'trouble' occurs. (At CLAS it is about 2 GeV 2 .) 

Being sure that one really deals with the leading twist contribution, such differences 
are not dramatic — as long as we are not interested in a high precision study of the effect. 
But how can we a priori be sure that there are no power corrections? Little is known 
about such corrections to the Sivers effect. For illustrative purposes, let us assume that 

A uf m JsLd = {'twist-2 Sivers effect' in Eqs. (11, 15) j +C(Q) (22) 

as one generically may expect. The 'coefficient' C(Q) could, in general, be flavour- 
dependent and typically depend on scale logarithmically (and depend on x, z, ... etc.). 
By looking at (21, 22) we see, that power corrections — if they play a role — are about 
5 times smaller at COMPASS compared to HERMES at x ~ 0. 15. 

Maybe such corrections are irrelevant for Q 2 > 1 GeV 2 which is typically used as DIS- 
cut. In any case, a careful comparison of all (present and future) data from COMPASS, 
HERMES and JLab will shed light on the possible size of power corrections. 

In particular, data from CLAS (first transverse target data will be available in 2011 
[86]) will be valuable — were a wide kinematical range is covered with very high 
statistics. Moreover, even after the 12 GeV beam energy upgrade, a certain fraction of the 
beam time also somehow lower beam energies will be available. Therefore, CLAS could 
provide precise information on SSAs at fixed value of x but for different Q 2 . Such data 
would allow to discriminate between 'leading twist' contributions and possible power 
corrections — paving the way to a concise understanding of the novel effects. 

One should not forget that in general 'power corrections' are by no means only 
'contaminations' to the parton model formalism. Rather, they also contain interesting 
physics, we refer to [87] for an example. 



11. CONCLUSIONS 



Since the release of first data on Sivers effect in SIDIS by HERMES and COMPASS 
[5, 6], more detailed, higher statistics (published and preliminary) data became available 
[7, 8, 9, 10, 11, 12, 13]. While the first data did not allow to constrain Sivers sea quarks 
[39, 40, 41, 42, 45], it is clear that the most recent ones [10, 11] can only be described, 
if the Sivers sea quarks are not zero. 

In this work we studied the recent preliminary data [10, 11]. For that we introduced 
relevant Sivers sea quarks (77, d, s, s), and choose a simple Ansatz that assumes the 
Sivers functions to be proportional to the unpolarized ones. The statistically satisfactory 
fit provides a good overall description of the preliminary data [10, 11]. 

The resulting Sivers u, d, 77, d distributions are in excellent agreement with large N c 
predictions [65]. The Sivers u, d distributions confirm earlier analyses [39, 40, 41, 42, 
45], and support the physical picture of the Sivers effect from [67]. The fit favors Sivers 
s, s distributions of opposite sign, and close to positivity bounds (in the model-dependent 
way we implemented them), though this observation has less statistical significance. 

That such a fit, that assumes f^j^ a (x) /f (x) for all flavours, works implies the 
following. The present data do not yet give much insights into details of the shapes 
of the Sivers functions. But they already tell us something about their magnitudes and 
relative signs, and this situation will improve due to the impact of future data. 

One could be happy about this situation. Earlier studies and the optimism, concerning 
measurability of the Sivers effect in DY, are confirmed [39, 40, 41, 42, 45]. The results 
fit into theoretical expectations [65, 67]. There is a rough, qualitative agreement with 
model results [88, 89, 90, 91, 92, 93] especially concerning Sivers quark distributions. 
But there is also a grain of salt. 

The description of the preliminary data [10] on the K + Sivers effect around x ~ 0.15 
is not ideal. Actually, we speak about 2 out of 61 data points 3 available presently on the 
Sivers effect in DIS production of various hadrons from various targets. Could these two 
data points be statistical fluctuations? Should they make us worry? Could they indicate 
power corrections? Or do they hint at novel effects? 

At the present stage it is not possible to draw definite conclusions. Unfortunately, 
the HERMES experiment is terminated. But fortunately there are promising prospects 
due to COMPASS with a proton target, and JLab where different transversely and 
longitudinally polarized targets will be explored [86]. 

The data from these experiments, which are complementary from the point of view 
of kinematics, will help to answer these and many other questions — on the Sivers and 
many other effects [41, 94, 95, 96, 97, 92, 98]. 



i We mean here data on x-dependence from [10, 11] (data on neutral kaons from [13] are not included in 
this counting). The relevant events, when binned correspondingly, also 'pop up' in data on, for example, 
the z-dependence and cause 'mismatches' also there, see Sec. 7. But these are not 'additional mismatches', 
since the different data sets are correlated. 
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